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(54) Lithographic apparatus 

(57) A device manufacturing method is disclosed in 
which the aberration of the projection system of a litho- 
graphic projection apparatus is obtained in terms of the 
Zernike expansion. The field distribution of displace- 
ment error and focal plane distortion of the projected im- 
age are calculated on the basis of the Zernike aberration 
and sensitivity coefficients which quantify the relation- 
ship between Zernike aberration components and the 
error in the image. A calculation is then performed to 
determine the compensation to apply to the apparatus 
in order to minimize the error in the image. The compen- 
sation is then applied to the apparatus. The compensa- 
tion may comprise increasing one component of aber- 
ration of the apparatus in order to decrease the effect of 
another aberration, such that, on balance, the image 
quality as a whole is improved. 
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Description 

[0001] The present invention relates to a lithographic projection apparatus comprising: 

5 - a radiation system for supplying a projection beam of radiation; 

a support structure for supporting patterning means, the patterning means serving to pattern the projection beam 

according to a desired pattern; 

a substrate table for holding a substrate; and 

a projection system for projecting the patterned beam onto a target portion of the substrate. 

10 

[0002] The term "patterning means" as here employed should be broadly interpreted as referring to means that can 
be used to endow an incoming radiation beam with a patterned cross-section, corresponding to a pattern that is to be 
created in a target portion of the substrate; the term "light valve" can also be used in this context. Generally, the said 
pattern will correspond to a particular functional layer in a device being created in the target portion, such as an inte- 
. 15 grated circuit or other device (see below). Examples of such patterning means include: 

A mask/The concept of a mask is well known in lithography, and it includes mask types such as binary, alternating 
phase-shift, and attenuated phase-shift, as well as various hybrid mask types. Placement of such a mask in the 
radiation beam causes selective transmission (in the case of a Iransmissive mask) or reflection (in the case of a 

20 reflective mask) of the radiation impinging on the mask, according to the pattern on the mask. In the case of a 

mask, the support structure will generally be a mask table, which ensures that the mask can be held at a desired 
position in the incoming radiation beam, and that it can be moved relative to the beam if so desired. 
A programmable mirror array. One example of such a device is a matrix-addressable surface having a viscoelastic 
control layer and a reflective surface. The basic principle behind such an apparatus is that (for example) addressed 

25 areas of the reflective surface reflect incident light as diffracted light, whereas unaddressed areas reflect incident 

light as undiffracted light. Using an appropriate filter, the said undiffracted light can be filtered out of the reflected 
beam, leaving only the diffracted light behind; in this manner, the beam becomes patterned according to the ad- 
dressing pattern of the matrix-adressable surface. An alternative embodiment of a programmable mirror array 
employs a matrix arrangement of tiny mirrors, each of which can be individually tilted about an axis by applying a 

30 suitable localized electric field, or by employing piezoelectric actuation means. Once again, the mirrors are matrix- 

addressable, such that addressed mirrors will reflect an incoming radiation beam in a different direction to unad- 
dressed mirrors; in this manner, the reflected beam is patterned according to the addressing pattern of the matrix- 
adressable mirrors. The required matrix addressing can be performed using suitable electronic means. In both of 
the situations described hereabove, the patterning means can comprise one or more programmable mirror arrays. 

35 More information on mirror arrays as here referred to can be gleaned, for example, from United States Patents 

US 5,296,891 and US 5,523,193, and PCT patent applications WO 98/38597 and WO 98/33096, which are incor- 
porated herein by reference. In the case of a programmable mirror array, the said support structure may be em- 
bodied as a frame or table, for example, which may be fixed or movable as required. 

A programmable LCD array. An example of such a construction is given in United States Patent US 5,229,872, 
40 which is incorporated herein by reference. As above, the support structure in this case may be embodied as a 

frame or table, for example, which may be fixed or movable as required. 

For purposes of simplicity, the rest of this text may, at certain locations, specifically direct itself to examples involving 
a mask and mask table; however, the general principles discussed in such instances should be seen in the broader 

45 context of the patterning means as hereabove set forth. 

[0003] Lithographic projection apparatus can be used, for example, in the manufacture of integrated circuits (ICs). 
In such a case, the patterning means may generate a circuit pattern corresponding to an individual layer of the IC. and 
this pattern can be imaged onto a target portion (e.g. comprising one or more dies) on a substrate (silicon wafer) that 
has been coated with a layer of radiation-sensitive material (resist). In general, a single wafer will contain a whole 

so network of adjacent target portions that are successively irradiated via the projection system, one at a time. In current 
apparatus, employing patterning by a mask on a mask table, a distinction can be made between two different types of 
machine. In one type of lithographic projection apparatus, each target portion is irradiated by exposing the entire mask 
pattern onto the target portion in one go; such an apparatus is commonly referred to as a wafer stepper. In an alternative 
apparatus — commonly referred to as a step-and-scan apparatus — each target portion is irradiated by progressively 

55 scanning the mask pattern under the projection beam in a given reference direction (the "scanning" direction) while 
synchronously scanning the substrate table parallel or anti-parallel to this direction; since, in general, the projection 
system will have a magnification factor M (generally < 1), the speed V at which the substrate table is scanned will be 
a factor M times that at which the mask table is scanned. More information with regard to lithographic devices as here 
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described can be gleaned, for example, from US 6,046,792, incorporated herein by reference. 
[0004] In a manufacturing process using a lithographic projection apparatus, a pattern (e.g. in a mask) is imaged 
onto a substrate that is at least partially covered by a layer of radiation-sensitive material (resist). Prior to this imaging 
step, the substrate may undergo various procedures, such as priming, resist coating and a soft bake. After exposure, 

5 the substrate may be subjected to other procedures, such as a post-exposure bake (PEB), development, a hard bake 
and measurement/inspection of the imaged features. This array of procedures is used as a basis to pattern an individual 
layer of a device, e.g. an IC. Such a patterned layer may then undergo various processes such as etching, ion-implan- 
tation (doping), metallization, oxidation, chemo-mechanical polishing, etc., all intended to finish off an individual layer. 
If several layers are required, then the whole procedure, or a variant thereof, will have to be repeated for each new 

10 layer. Eventually, an array of devices will be present on the substrate (wafer). These devices are then separated from 
one another by a technique such as dicing or sawing, whence the individual devices can be mounted on a carrier, 
connected to pins, etc. Further information regarding such processes can be obtained, for example, from the book 
"Microchip Fabrication: A Practical Guide to Semiconductor Processing", Third Edition, by Peter van Zant, McGraw 
Hill Publishing Co., 1997, ISBN 0-07-067250-4, incorporated herein by reference. 

js [0005] For the sake of simplicity, the projection system may hereinafter be referred to as the "lens"; however this 
term should be broadly interpreted as encompassing various types of projection system, including refractive optics, 
reflective optics, and catadioptric systems, for example. The radiation system may also include components operating 
according to any of these design types for directing, shaping or controlling the projection beam of radiation, and such 
components may also be referred to below, collectively or singularly, as a "lens". Further, the lithographic apparatus 

20 may be of a type having two or more substrate tables (and/or two or more mask tables). In such "multiple stage" devices 
the additional tables may be used in parallel, or preparatory steps may be carried out on one or more tables while one 
or more other tables are being used for exposures. Dual stage lithographic apparatus are described, for example, in 
US 5,969,441 and WO 98/40791 , incorporated herein by reference. 

[0006] When performing imaging in a lithographic projection apparatus, despite the groat care with which the pro- 

25 jection system is designed and the very high accuracy with which the system is manufactured and controlled during 
operation, the image can still be subject to aberrations such as, for example, distortion (i.e. a non-uniform image 
displacement in the target portion at the image plane: the XY-plane), lateral image shift (i.e. a uniform image displace- 
ment in the target portion at the image plane), image rotation, asymmetric magnification, and focal plane deformation 
(i.e. a non-uniform image displacement in the Z-direction, for instance due to field curvature). Notice that, in general, 

30 image errors are not necessarily uniform, and can vary as a function of position in the image field. Distortion and focal 
plane deformation can lead to overlay and focus errors, for example overlay errors between different mask structures, 
and line-width errors. As the size of features to be imaged decreases, these errors can become intolerable. 
[0007] Consequently, it is desirable to provide compensation (such as adjustment of the projection system and/or 
substrate) to correct for, or at least attempt to minimize, these errors. This presents the problems of first measuring 

35 the errors and then calculating appropriate compensation. Previously, alignment systems were used to measure the 
displacements in the image field of alignment marks. However, alignment marks typically consist of relatively large 
features (of the order of a few microns), causing them to be very sensitive to aberrations of the projection system. The 
alignment marks are unrepresentative of the actual features being imaged, and because the imaging errors depend 
inter alia on feature size, the displacements measured and compensations calculated did not necessarily optimize the 

40 image for the desired features. 

[0008] Another problem occurs when, for instance because of residual manufacturing errors, the projection system 
features an asymmetric variation of aberration over the field. These variations may be such that at the edge of the field 
the aberration becomes intolerable. 

[0009] A further problem occurs when using phase-shift masks (PSM's). Conventionally, the phase shift in such 
45 masks has to be precisely 1 80 degrees. The control of the phase is critical; deviation from 180 degrees is detrimental. 
PSM's, which are expensive to make, must be carefully inspected, and any masks with substantial deviation in phase 
shift from 180 degrees will generally be rejected. This leads to increased mask prices. 

[001 0] A further problem occurs with the increasing requ irements imposed on the control of critical dimension ("CD"). 
The critical dimension is the smallest width of a line or the smallest space between two lines permitted in the fabrication 

so of a device. In particularthe control of the uniformity of CD, the so-called "CD uniformity", is of importance. In lithography, 
efforts to achieve better line width control and CD uniformity have recently led to the definition and study of particular 
error types occurring in features, as obtained upon exposure and processing (see description above). For instance, 
such image error types are an asymmetric distribution of CD over a target portion, an asymmetry of CD with respect 
to defocus (which results in a tilt of Bossung curves), asymmetries of CD within a feature comprising a plurality of bars 

55 (commonly referred to as Left-Right asymmetry), asymmetries of CD within a feature comprising either two or five bars 
(commonly known as L1-L2 and L1 -L5, respectively), differences of CD between patterns that are substantially directed 
along two mutually orthogonal directions (for instance the so-called "H-V" lithographic error), and for instance a variation 
of CD within a feature, along a bar, commonly known as "C-D". Just as the aberrations mentioned above, these errors 
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are generally non-uniform over the field. For simplicity we will hereafter refer to any of these error types - including the 
errors such as, for example, distortion, lateral image shift, image rotation, asymmetric magnification, and focal plane 
deformation - as "lithographic errors", i.e. feature-deficiencies of relevance for the lithographer. 
[0011] Lithographic errors are caused by specific properties of the lithographic projection apparatus. For instance, 

5 the aberration of the projection system, or imperfections of the patterning means and imperfections of patterns gener- 
ated by the patterning means, or imperfections of the projection beam may cause lithographic errors. However, also 
nominal properties (i.e. properties as designed) of the lithographic projection apparatus may cause unwanted litho- 
graphic errors. For instance, residual lens aberrations which are part of the nominal design may cause lithographic 
errors. For reference hereafter, we will refer to any such properties that may cause lithographic errors as "properties". 

10 [001 2] As mentioned above } the image of a pattern can be subject to aberrations of the projection system. A resulting 
variation of CD (for example, within a target portion) can be measured and subsequently be mapped to an effective 
aberration condition of the projection system which could produce said measured CD variation. A compensation can 
than be provided to the lithographic projection system such as to improve CD uniformity. Such a CD-control method 
is described in United States Patent US 6,1 15,1 08, incorporated herein by reference, and comprises imaging a plurality 

is of test patterns at each field point of a plurality of field points, a subsequent processing of the exposed substrate, and 
a subsequent CD measurement for each of the imaged and processed test patterns. Consequently, the method is time 
consuming and not suitable for in-sltu CD control. With increasing demands on throughput {i.e. the number of substrates 
that can be processed in a unit of time) as well as CD uniformity, the control, compensation and balancing of lithographic 
errors must be improved, and hence, there is the problem of furthering appropriate control of properties. 

20 [0013] It is an object of the present invention to alleviate the above problems. 

[001 4] This and other objects are achieved according to the invention in a device manufacturing method comprising 
the steps of: 

providing a substrate that is at least partially covered by a layer of radiation-sensitive material; 
providing a projection beam of radiation using a radiation system; 

using patterning means to endow the projection beam with a pattern in its cross-section; 

projecting the patterned beam of radiation onto a target portion of the layer of radiation-sensitive material using a 
projection system, 

30 characterized by: 

obtaining information on properties of at least one of the substrate, the layer of radiation sensitive material, the 
projection beam, the patterning means and the projection system; 

obtaining a plurality of coefficients which quantify the relationship between said properties and at least one of a 
35 plurality of lithographic errors causing anomaly in a projected image in the radiation sensitive layer; 

defining a merit function which weighs and sums lithographic errors; 

calculating a compensation to apply to at least one of the substrate, the projection beam, the patterning means 
and the projection system to optimize the merit function; and 
applying the calculated compensation. 

40 

[0015] Although preferably all lithographic errors shall be minimized, this will in general not be possible since there 
are not enough compensators available. Therefore, it is useful to define a merit function which suitably describes the 
quality of the imaging process in terms of a sum of weighted lithographic errors. A weight assigned to a lithographic 
error shall be indicative of the relative importance of that lithographic error. Compensation can be used to optimize (i. 

45 e . to minimize the value of) the merit function, such as to optimize the imaging and the resulting lithographic process. 
As explained above : lithographic errors depend on properties (of the lithographic projection apparatus) and can be 
calculated using, for instance, commercially available lithography simulation software such as Prolilh™, Solid-C™ or 
LithoCruiser ™. For instance, given specific (critical) pattern features to be imaged, given specific pattern errors (such 
as phase errors with a phase shift mask), given the aberration of the projection system, given the data concerning the 

50 radiation sensitive layer on the substrate, and given the radiation beam properties such as radiation energy and wave- 
length (i.e. given a set of properties), predictions regarding CD uniformity or feature-specific deficiencies can be made 
with these simulation programs. By introducing small variations of the properties (for instance, by introducing small 
variations of the projection lens aberration) and calculating the corresponding changes in the lithographic errors, co- 
efficients quantifying the relationship between a lithographic error and said properties can be established. Obtaining 

55 information on properties comprises (in-situ) measuring of properties. For example, the actual aberration condition of 
the projection system can be measured in-situ using a suitable technique such as one of those disclosed in European 
Patent Application Number EP 1 128 217 A2 and in P. Venkataraman, et al., "Aberrations of steppers using Phase 
Shifting Point Diffraction lnterferometry ,, , in Optical MicroUthography XII t, J. Progler, Editor, Proceedings of SPIE Vol. 
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4000, 1 245-1249 (2000). An advantage of the method according to the present invention is that the method does not 
necessarily comprise the step of processing an exposed substrate. 

[001 6] Key to the control of lithographic errors is the ability to compensate, or at least affect, these lithographic errors. 
Any means enabling such control will be referred to hereafter as compensation means. By the term "compensation" 

5 mentioned above any change of a property caused by activating a compensation means is meant, and we will use the 
term compensation in this sense hereafter. Compensation means suitable for use with a lithographic projection appa- 
ratus are, for instance, means enabling a fine positioning (an X-, Y-, and Z- translation, and a rotation about the Y-, 
and Z-axis) of the holder for holding the patterning means, means enabling a similar fine positioning of the substrate 
table, means to move or deform optical elements (in particular, to fine position, using an X-, Y-, and Z- translation/ 

10 rotation, optical elements of the projection system), and, for instance, means to change the energy of the radiation 
impinging on the target portion. However, suitable compensation means are not limited to said examples: for instance, 
means to change the wavelength of the radiation beam, means to change the pattern, means to change the index of 
refraction of gas-filled spaces traversed by the projection beam, and means to change the spatial distribution of the 
intensity of the radiation beam may also serve to affect lithographic errors. 
.15 [0017] The optimization calculation is based on varying the available compensations. The calculation comprises 
minimization of the merit function. Since the merit function is available as a mathematical expression (relating litho- 
graphic errors to compensations), standard minimization routines such as a "weighted least squares optimization" can 
be used to find a set of compensations at which the merit function has a minimum. 

[0018] By performing the optimization calculation, the quality of the imaging can be objectively improved, as defined 
20 by the merit function. The compensation enables, for example, one component of aberration of the projection system 
to be increased in order to decrease the effect of another aberration component, such that, on balance, the image 
quality as a whole is improved. In other words, preferably the improvement in one aspect of the imaging more than 
compensates for a change of a different aspect of the imaging. According to another example, the patterning means 
can bo translated, rotated and/or tilted (a form of rotation) in order to minimize the overall local displacement of the 
25 image from it's ideal position across the whole image field. In a further example, linear comatic aberration can be 
introduced within the projection system to reduce or eliminate left-right asymmetry caused by 3-wave aberration. In a 
further example, spherical aberration can be introduced to compensate for "Bossung tilt" caused by a phase shift error 
in a phase shift mask, which would otherwise have to be rejected. Bossung tilt is the lithographic error whereby at best 
focus position there is a gradient in the graph of exposed feature width against focus position (said graph being a so- 
30 called Bossung curve). 

[0019] Said properties (of, for example, the projection system) can be stored in a database for use in optimization 
calculations for exposures using the apparatus containing that projection system. The coefficients which quantify the 
relationship between those properties and a lithographic effect can also be stored in a database as sets or families of 
coefficients depending on pattern feature type, size, orientation, illumination mode, numerical aperture and so on. 
35 [0020] According to a further aspect of the invention there is provided a lithographic projection apparatus comprising: 

a radiation system for providing a projection beam of radiation; 

a support structure for supporting patterning means, the patterning means serving to pattern the projection beam 
according to a desired pattern; 
40 - a substrate table for holding a substrate; 

a projection system for projecting the patterned beam onto a target portion of the substrate, 

characterized by further comprising: 

45 compensation means appliable to at least one of the holder for holding the patterning means, the substrate table, 

the radiation system, the patterning means and the projection system, 

to optimize a merit function which weighs and sums lithographic errors causing anomaly in a projected image in 
the radiation sensitive layer; and 

a processor for calculating at least one compensation to be applied by said compensation means on the basis of 
so a plurality of coefficients which quantify the relationship between at least one lithographic error and properties of 

at least one of the patterning means, the projection system, the radiation sensitive layer on the substrate and the 
projection beam. 

[0021] Although specific reference may be made in this text to the use of the apparatus according to the invention 
55 in the manufacture of ICs, it should be explicitly understood that such an apparatus has many other possible applica- 
tions. For example, it may be employed in the manufacture of integrated optical systems, guidance and detection 
patterns for magnetic domain memories, liquid-crystal display panels, thin-film magnetic heads, etc. The skilled artisan 
will appreciate that, in the context of such alternative applications, any use of the terms "reticle 0 , "wafer" or "die" in this 
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text should be considered as being replaced by the more general terms "mask", "substrate" and "target portion 0 , re- 
spectively. 

[0022] In the present document, the terms "radiation" and "beam" are used to encompass all types of electromagnetic 
radiation, including ultraviolet radiation {e.g. with a wavelength of 365, 248, 193, 157 or 126 nm) and EUV (extreme 
5 ultra-violet radiation, e.g. having a wavelength in the range 5-20 nm). 

[0023] Embodiments of the invention will now be described, by way of example only, with reference to the accom- 
panying schematic drawings in which: 

Figure 1 depicts a lithographic projection apparatus according to an embodiment of the invention; and 
10 Figure 2 depicts the imaging of a decentered pattern, where the decenter is the compensator for alleviating the 

effects of a lithographic error. The graph at the bottom shows the boundary of the area of x, yf ield points contributing 
to the merit function, in relation to the boundary of the full field target portion. The horizontal axis represents the 
x-position along the slit, the vertical axis represents the y-position (a position along the scanning direction). The 
graph in the middle shows a plot of a lithographic error, which is asymmetrically distributed along the x-axis. The 
15 horizontal axis represents the x-position along the slit, the vertical axis represents the magnitude of a lithographic 

error. 

Embodiment 1 

20 [0024] Figure 1 schematically depicts a lithographic projection apparatus according to a particular embodiment of 
the invention. The apparatus comprises: 

a radiation system Ex, IL, for supplying a projection beam PB of radiation (e.g. UV radiation such as for example 
generated by an excimer laser operating at a wavelength of 24B nm, 1 93 nm or 1 57 nm, or by a laser-fired plasma 
25 source operating at 13,6 nm). In this particular case, the radiation system also comprises a radiation source LA; 

a first object table (mask table) MT provided with a mask holder for holding a mask MA (e.g. a reticle), and connected 
to first positioning means for accurately positioning the mask with respect to item PL; 

a second object table (substrate table) WT provided with a substrate holder for holding a substrate W (e.g. a resist- 
coated silicon wafer), and connected to second positioning means for accurately positioning the substrate with 
30 respect to item PL; 

a projection system ("lens") PL (e.g. a quartz and/or CaF 2 lens system or a catadioptric system comprising lens 
elements made from such materails, or a mirror system) for imaging an irradiated portion of the mask MA onto a 
target portion C (e.g. comprising one or more dies) of the substrate W. 

35 As here depicted : the apparatus is of a transmissive type (i.e. has a transmissive mask). However, in general, it may 
also be of a reflective type, for example (with a reflective mask). Alternatively, the apparatus may employ another kind 
of patterning means, such as a programmable mirror array of a type as referred to above. 

[0025] The source LA (e.g. a UV excimer laser, a laser-fired plasma source, a discharge source, or an undulator or 
wiggler provided around the path of an electron beam in a storage ring or synchrotron) produces a beam of radiation. 

40 This beam is fed into an illumination system (illuminator) IL, either directly or after having traversed conditioning means, 
such as a beam expander Ex, for example. The illuminator IL may comprise adjusting means AM for setting the outer 
and/or inner radial extent (commonly referred to as o-outer and o-inner, respectively) of the intensity distribution in the 
beam. In addition, it will generally comprise various other components, such as an integrator IN and a condenser CO. 
In this way, the beam PB impinging on the mask MA has a desired uniformity and intensity distribution in its cross- 

45 section. 

[0026] It should be noted with regard to Figure 1 that the source LA may be within the housing of the lithographic 
projection apparatus (as is often the case when the source LA is a mercury lamp, for example), but that it may also be 
remote from the lithographic projection apparatus, the radiation beam which it produces being led into the apparatus 
(e.g. with the aid of suitable directing mirrors); this latter scenario is often the case when the source LA is an excimer 

50 laser. The current invention and Claims encompass both of these scenarios. 

[0027] The beam PB subsequently intercepts the mask MA, which is held on a mask table MT. Having traversed the 
mask MA, the beam PB passes through the lens PL, which focuses the beam PB onto a target portion C of the substrate 
W. With the aid of the second positioning means (and interferometric measuring means IF), the substrate table WT 
can be moved accurately, e.g. so as to position different target portions C in the path of the beam PB. Similarly, the 

55 first positioning means can be used to accurately position the mask MA with respect to the path of the beam PB, e.g. 
after mechanical retrieval of the mask MA from a mask library, or during a scan. In general, movement of the object 
tables MT, WT will be realized with the aid of a long-stroke module (course positioning) and a short-stroke module (fine 
positioning), which are not explicitly depicted in Figure 1. However, in the case of a wafer stepper (as opposed to a 
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step-and-scan apparatus) the mask table MT may just be connected to a short stroke actuator, or may be fixed. 
[0028] The depicted apparatus can be used in two different modes: 

1. In step mode, the mask table MT is kept essentially stationary, and an entire mask image is projected in one go 
5 (i.e. a single "flash") onto a target portion C. The substrate table WT is then shifted in the x and/or y directions so 

that a different target portion C can be irradiated by the beam PB; 

2. In scan mode, essentially the same scenario applies, except that a given target portion C is not exposed in a 
single "flash", instead, the mask table MT is movable in a given direction (the so-called "scan direction", e.g. the 
y direction) with a speed v, so that the projection beam PB is caused to scan over a mask image; concurrently, 

10 the substrate table WT is simultaneously moved in the same or opposite direction at a speed V= Mv, in which M 

is the magnification of the lens PL (typically, M= 1/4 or 1/5). In this manner, a relatively large target portion C can 
be exposed, without having to compromise on resolution. 

[0029] The image distortion and focal plane deformation occurring in the projected image depend on various factors, 
is such as the lens aberration, illumination mode used [e.g. numerical aperture (NA), pupil filling (sigma setting)), and the 
type and size of the feature imaged. The lens aberration can be expressed in terms of the Zernike expansion, in which 
the aberrated wavefront is given by the sum of a plurality of Zemike polynomials (functions of position in the image 
field) each multiplied by a respective Zernike coefficient. The image X, Y and Z displacements can be modeled by the 
following linear expressions: 

20 

dX{x.y)^^Zj{x,y) (1) 



dY{x,y) = ^~Zj(x,y) (2) 



30 j dZj 



35 



40 where: 



dF{x,y) = Y^' Z j^>y) < 3 > 



dX, dY, dF are the respective X, Y and Z displacements; 
(x, y) is a specific field location; 
j is an index for the jth Zernike aberration; 
45 Zj is an aberration coefficient representative of the contribution of the yth Zernike polynomial to the wavefront 

aberration at field point (x,y); and 

5— is the partial derivative of the displacement with respect to the jth Zernike aberration. 

[0030] Notice that dX, dVand dFare lithographic errors, and that the coefficients Z y - are properties (of the imaging 
so system including an object and image plane). 

[0031] A further refinement is to write the Z-displacement dFas two components oVand dH which are the optimal 
focus positions along the Z-axis for so-called vertical and horizontal lines, respectively; these are features extending 
along the yand xaxes, respectively. Using this, and a different notation, the above equations can be rewritten, for the 
37 Zernike aberrations, as: 



55 



dX(x,y) = ,41 ,Z1 (x,y) + A2.Z2(x f y) + ... + A37.Z37{x,y) (4) 
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dY(x,y) = B1.Z1(x,y) + B2.Z2(x,y) + ... + B37.Z37(x,y) (5) 
dH(x,y) = C1 ,Z1 (x,y) + C2.Z2(x,y) +...+ C37.Z37(x,y) (6) 
a^(x,y) = D1 .Z1 (x,y) + D2.Z2(x,y) +...+ D37.Z37(x,y) (7) 



to where the partial derivatives have been replaced by equivalent sensitivity coefficients Ai, Bi, Ci, Di which express inter 
alia the effects of illumination mode, NA, sigma and pattern feature size and type. 
[0032] The procedure according to this embodiment of the invention is as follows: 

1 . The intrinsic aberration of the lens is measured using a suitable technique, such as one of those disclosed in 
15 European Patent Application EP 1 128 217 A2 ; to obtain the Zernike coefficients (or previously measured values 

can be used). The Zernike coefficients are then multiplied by their respective Zernike functions to obtain a "field 
map" of the Zernike aberrations (Z/(x,y) in the above notation), which is stored. This information depends only on 
the lens and is independent of the illumination and features. 

2. For the exposure of a specific pattern feature, at a particular illumination setting, the sensitivity coefficients Ai, 
20 Bi, Ci, Di are calculated using known computational techniques, such as image simulation software, for exampie 

the program known as So!id-C™ supplied by Sigma-C GmbH of Germany. Other suitable software packages, such 
as one known as Proiith™ or LithoCruiser™, may alternatively be used. These sensitivity coefficients do not depend 
on the field location, so only need to be calculated once per feature and illumination setting combination, and are 
then stored in a database. 

25 3. For the exposure of a specific pattern feature, at a particular illumination setting, the field distribution of the X, 

Y and Z-displacements is calculated by means of the models, for example equations (4) to (7), and using the field 
map from step 1 and the sensitivity coefficients from step 2. A calculation is performed to determine what com- 
pensation to provide to minimize the image displacements for the whole field. Such a calculation can be based on 
a merit function of the form 

30 

S = W iZ(dXix,y)Y + wl^idY^y)) 2 + w3^(dH (x, y)) 2 + ^{dV^y)) 1 (8) 

x % y x,y x,y x,y 

35 

where wi(i= 1,2,3,4) are weights, and the sum over an appropriately chosen grid of field coordinates x,y ensures an 
averaging of the lithographic error over a target portion. The weights can be chosen at the discretion of the lithographer. 
For dX, dY, dH and oVone can substitute expressions (4), (5), (6) and (7) respectively. A simple compensator is, for 
40 instance, the Z-position of the substrate, commonly called "focus". It is well known that a change of focus, by moving 
the substrate through a distance dZ, affects primarily the value of Z4, and in this particular case, the error is independent 
of the x,y position in the field. So, in order to incorporate the error of the compensation dZon the merit function (through 
the errors that dZ has on dX, dY, dH and dV), one has to substitute Z4 + F4(dZ,x,y) for Z4 in the expressions (4), (5), 
(6) and (7), as used in equation (8), where the function F4(dZ,x,y) is simply given by 

45 

F4(dZ,x,y) = a4-dZ (9) 

a4 being a (field independent) constant of proportionality which follows from aberration theory. Notice, however, that 
so in practice the compensation dZwill, besides this primary impact on Z4, also affect, as a secondary error, other Zernike 
aberration coefficients. Generally, the error of a compensation C comprises a primary error on a few Zernike coefficients 
and a secondary change of balance between the other Zernike aberration coefficients. So, the error of a set of n 
compensations Ci, with / =1 ,2,..,n, on the merit function can be visualized by replacing Z/(x,y) in the expressions (4), 
(5), (6) and (7), as used in equation (8), by Z/(x,y,C1 ,C2 C\,...,Cn), where 

55 
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Zj(x 9 >•, CI, C2, • ■ , Ci, • - • , Of) = Zy(;c, y) + £ F/(Cz, x, >) . (10) 



to 



f5 



20 



SO 



1=1 



The functions Fj(Ci,x,y) are (to good approximation) linearly dependent on C/'and vanish for C*=0. Therefore, to a good 
approximation one obtains: 



Zj{x 9 y,C\ 9 C2 n •-sCi^--,C«).*^.>') + £CiLyi(*.y) (H) 



The constants of proportionality Uy7(x,y) are defined by 

Uii(x ,y )= mCLm, (12) 



and follow from the theory of aberrations or can be calculated using optical design software, such as for instance 
CodoV™. Notice that, in general, the constants of proportionality are field dependent. In the example above, a4 rep- 

25 resents a special case where the constant of proportionality is field independent. 

[0033] In view of the above, the merit function S depends on the compensations CI, C2 : .... Cn. Then, using a 
standard minimization procedure such as a weighted least square optimization (see for instance the publication "Nu- 
merical Recipes in C", First Edition, by William H. Press, Brian P. Flannery, Saul A. Teukolsky, and William T. Vetterling, 
Cambridge University Press, 1988), values of the compensations can be found for which S is at a minimum value. 

30 Next, the projection apparatus is adjusted to provide the calculated compensation. Examples of available compensa- 
tions include: 

translation of the mask table (to shift the x, y co-ordinates); 

lens adjustment (e.g. to change the magnification and/or field curvature, using internal lens-element manipulators 
35 and/or telecentricity adjustment); 

rotation of the mask (reticle); 

adjustment of the mask height (e.g. to compensate for 3 rd order distortion); 
adjustment of the mask tilt; and 

adjustment of leveling (e.g. to compensate for z-displacement and tilt). 

40 

[0034] For a particular pattern to be imaged, the illumination setting to be used is usually predetermined to optimize 
contrast. Similarly, although the pattern may have many feature types, in practice one feature type will be the most 
critical. At step 2, the sensitivity coefficients can be calculated for the predetermined illumination setting and the most 
critical feature type. Alternatively, at step 2, a family of sets of sensitivity coefficients can be derived for a range of 
45 feature types and illumination settings and stored in a database. Then the set of sensitivity coefficients for a pattern, 
with particular most critical feature type and predetermined illumination setting, can simply be looked up in the database 
of coefficient families. The sensitivity coefficients, derived either way, can then be stored as the default values for that 
particular pattern. If the model of equations (4) to (7) is used there will only be 37x4 coefficient values associated with 
each pattern. 



Embodiment 2 

[0035] As described above, a translation of the mask table (holder for the patterning means) is one of the compen- 
sations available. In one embodiment this translational degree of freedom is exploited to establish a decentered imaging 
of a pattern: see Fig. 2. A projection system PL of a scanning lithography apparatus is shown in Fig. 2 and is assumed 
to exhibit a strong asymmetric distribution 21 of, for instance, Z7(x,y) along the x-direction. Such an asymmetry may 
be caused by residual manufacturing errors that occurred during the manufacturing of the projection system, and is 
known to be an important contributor to the lithographic error Left-Right asymmetry (LR asymmetry). If the pattern 23 
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that is to be imaged is smaller that the full field width 25 s one can shift the pattern in the x-direction over a distance 
231 , such that the patterned beam traverses the projection system in a decentered manner. As a result, the portion 27 
of the field which suffers strongly from Z7 is avoided, and the lithographic error LR asymmetry will be reduced. In this 
example, the merit function S will comprise a weighted term which accounts for LR asymmetry (with a weight wchosen 
by the lithographer): 



5 = ... 4- w^T {LRasymmetrytx, yjf + • • ■ (13) 



The compensation 231 affects, in this embodiment, the relative position of the boundary 29 of the grid of x,y field points 
(that contribute to the merit function) with respect to the boundary 291 of the full field target portion. 



[0036] As described in Embodiment 1 above, displacement-compensation means are used to minimize the X, Y and 
Z image displacements resulting from, for instance, aberration in the lens. Effectively the plane of the image and the 
wafer is shifted to balance aberration effects of the lens. Alternatively, or in addition, the aberration of the lens itself 
can be adjusted. For example, spherical aberration adjustment and comatic aberration adjustment can be provided 
(spherical and comatic aberration are the iowest order even and odd aberrations, respectively). In other respects, this 
embodiment is the same as Embodiment 1 . 

[0037] Aberration adjustment will influence the distribution of all aberrations (i.e. change the values of the Zernike 
coefficients), so will need to be taken into account at step 1 . The means of adjusting the aberration may also influence 
the sensitivity coefficients Ai, Bi, Ci, Di, so may need to be taken into account in step 2. Lens manipulation primarily 
influences the low order Zernike aberrations (Z2, Z3, Z4, Z5). The impact of these aberrations on the imaging is dis- 
placement and defocus (astigmatism), and this impact is independent of the pattern feature structure and the illumi- 
nation settings. 

[0038] Often it is best simply to minimize the aberration of the lens. However, as can be seen from equations (4) to 
(7), the ability to adjust the aberrations Zj, gives more adjustable parameters, and hence more degrees of freedom, 
for minimizing the displacement errors. 



[0039] The invention is not limited to reducing overlay and focus errors (minimizing the displacements). Compensa- 
tion can be used to reduce any lithographic effect that is detrimental to image quality. For example left-right asymmetry 
(LR asymmetry) caused by 3-wave aberration can be reduced by introducing linear comatic aberration within the sys- 
tem. In other words, increasing one aberration can be used to decrease another aberration such that, on balance, the 
quality as a whole of a printed feature is improved. It is then necessary to include LR asymmetry in the merit function. 
Let us denote an arbitrary lithographic error, that is to be minimized, by LEj,j= 1,2,3,4,.,. For instance, we may rewrite 
Equation (8) as 



if dX = LE1 , dY = LE2, dH = LE3 and dV = LE4. One can now include LR asymmetry in the merit function with a weight 
w5 t to be chosen by the lithographer: 




Embodiment 3 



Embodiment 4 



4 




(14) 



At \ 

S = Y, wjY.iLEj) 2 +w5^(LRasymmetty(x,y)) , 



(15) 



/=] \ *.y J x.y 
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With simulation software one can calculate the sensitivity coefficients for the lithographic error LRasymmetry(x,^ in 
analogy to equations (4) - (7), and proceed as in embodiment 1 . Another example of a lithographic error which can be 
compensated for is the error known as Bossung tilt (also known as "iso-focal tilt"). It may occur when a phase shift 
mask (PSM) is used as patterning means. Ideally, the phase shift between the regions of the mask is precisely 180 

5 degrees. Any deviation from 180 degrees phase shift will cause a Bossung tilt (i.e. even at the best focus position, 
there is a gradient in the graph of exposed feature width against focus position (Bossung curve); therefore the process 
latitude is very poor, and in practice one may have to discard the mask). However, it has been found that: (i) the phase 
error of a PSM; and(ii) spherical aberration (characterized by the Z9 Zernike parameter), have a similar impact on the 
Bossung curves, both leading to Bossung tilt. Therefore, a specific amount of spherical aberration can be introduced 

10 to cause Bossung tilt in one direction to cancel out a PSM phase error which causes Bossung tilt in the opposite 
direction. In this example, the merit function includes Bossung tilt, referred to hereafter as er(x,y), which is a function 
of lens aberrations; in analogy with, for instance, equation* (4) we may write 



15 



25 



35 



45 



BT(x,y) = P^ Z1 (x,y) + P2Z2(x,y) + - + P9Z9(x,y) + »• P37Z37(x,y) + O.G (16) 



In this example, the sensitivity coefficient P9 is relatively large with respect to the other sensitivity coefficients Py>1 ,2,.. 
8,10,. .,37, e is the PSM phase error, and Ois a sensitivity coefficient 

[0040] A compensation means which is known primarily to affect 29 can now be used for counteracting the litho- 
20 graphic error Bossung tilt. If an optimal compromise between correcting for Bossung tilt, LR asymmetry, and the image 
shifts mentioned in embodiment 1 is to be found, the merit function shall be chosen to be 



S = jTvJ^LEj) 2 ^wS^iLRasymmetryix^y)) 2 + v^X(£7X*,.y)) 2 (17) 



30 [0041] Whilst specific embodiments of the invention have been described above, it will be appreciated that the in- 
vention may be practiced otherwise than as described. The description is not intended to limit the invention. 



Claims 

1 . A device manufacturing method comprising the steps of: 



- providing a substrate that is at least partially covered by a layer of radiation-sensitive material; 
providing a projection beam of radiation using a radiation system; 

40 - using patterning means to endow the projection beam with a pattern in its cross-section; 

- projecting the patterned beam of radiation onto a target portion of the layer of radiation-sensitive material using 
a projection system, 



characterized by: 



obtaining information on properties of at least one of the substrate, the layer of radiation sensitive material, 
the projection beam, the patterning means and the projection system; 

obtaining a plurality of coefficients which quantify the relationship between said properties and at least one of 
a plurality of lithographic errors causing anomaly in a projected image in the radiation sensitive layer; 
so defining a merit function which weighs and sums lithographic errors; 

calculating a compensation to apply to at least one of the substrate, the projection beam, the patterning means 
and the projection system to optimize the merit function; and 
applying the calculated compensation. 

55 2. A device manufacturing method according to claim 1 , wherein said properties are selected from the group con- 
sisting of: pattern features; mask pattern errors; the aberration of the projection system; data concerning the ra- 
diation sensitive layer on the substrate; radiation beam energy and wavelength. 
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A device manufacturing method according to claim 1 or 2, wherein said at least one lithographic error is selected 
from the group consisting of: feature displacement in the X.Y-plane; feature rotation; asymmetric magnification of 
feature; distortion: focal plane deformation; asymmetric distribution of CD over a target portion; asymmetry of CD 
with respect to defocus; Left-Right asymmetries of CD; L1 -L2; L1 -L5; differences of CD between patterns that are 
substantially directed along two mutually orthogonal directions; variation of CD within a feature, along a bar, com- 
monly known as C-D; and combinations thereof. 

A device manufacturing method according to claim 1 , 2 or 3, wherein said coefficients comprise linear constants 
of proportionality between said properties and said lithographic errors. 

A device manufacturing method according to any one of the preceding claims, wherein said coefficients comprise 
a plurality of families of coefficients for at least one of: different illumination settings of the radiation system and/ 
or projection system; and different features of the pattern to be projected 

A device manufacturing method according to any one of the preceding claims, wherein the process of obtaining 
comprises at least one of: measuring; calculating; and reading from a database. 

A device manufacturing method according to any one of the preceding claims, wherein said properties comprise 
a field map of aberration as a function of position in the image to be formed at said target portion. 

A device manufacturing method according to any one of the preceding claims, wherein said calculation is performed 
to optimize the merit function for the most critical feature of the pattern to be projected. 

A device manufacturing method according to any one of the preceding claims, wherein the compensation com- 
prises at least one of: translation of the patterning means; rotation of the patterning means; translation of the 
substrate; rotation of the substrate; movement or deformation of optical elements, preferably optical elements of 
the projection system; change of energy of the radiation beam; change of wavelength of the radiation beam; change 
of index of refraction of gas-filled spaces traversed by the projection beam; change of the spatial distribution of 
the intensity of the radiation; change of used portion of field. 

10. A device manufacturing method according to any one of the preceding claims, wherein: 

said information on properties is a field map of aberration of the projection system expressed in terms of the 
Zernike expansion; 

35 said lithographic error comprises at least one of: the X displacement of the pattern; the Y displacement of the 

pattern; the Z displacement of features of the pattern oriented in the X direction; and the Z displacement of 
features of the pattern oriented in the Y direction, each as a function of position in the X, Y plane of the target 
portion of the substrate; and 

said coefficients are the constants of proportionality giving the sensitivity of the lithographic error to each 
40 Zernike aberration component in said expansion. 

11. A device manufacturing method according to any one of the preceding claims, wherein: 

said merit function describes the left-right asymmetry of the pattern; 
45 said coefficients relate the left-right asymmetry to the position of the pattern with respect to the full field of the 

projection system; and 

said compensation is adjustment of the position of the pattern with respect to the full field of the projection 
system. 

so 12. A device manufacturing method according claim 11 , wherein said compensation comprises a shift of the pattern 
in the X direction, such that the patterned beam traverses the projection system in a decentered manner. 

13. A device manufacturing method according to any one of the preceding claims, wherein: 

55 said merit function describes the left-right asymmetry of the pattern; 

said coefficients are the constants of proportionality giving the sensitivity of the left-right asymmetry to aber- 
ration of the projection system; and 

said compensation is adjustment of the linear comatic aberration of the projection system. 
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14. A device manufacturing method according to any one of the preceding claims, wherein: 
said patterning means comprises a phase shift mask; 

said information on properties comprises the phase shift error of the phase shift mask; 
5 said merit function describes Bossung tilt; 

said coefficients are the constants of proportionality giving the sensitivity of the Bossung tilt to the phase shift 

error and to aberration of the projection system; and 

said compensation is adjustment of the aberration of the projection system. 

10 15. A device manufacturing method according to claim 12, wherein said aberration of the projection system is partic- 
ularly the spherical aberration, and more preferably the Zernike aberration component known as Z9. 

16. A lithographic projection apparatus comprising: 

15 a radiation system for providing a projection beam of radiation; 

a support structure for supporting patterning means, the patterning means serving to pattern the projection 
beam according to a desired pattern; 
a substrate table for holding a substrate; 

a projection system for projecting the patterned beam onto a target portion of the substrate, 

20 

characterized by further comprising: 

compensation means appliable to at least one of the holder for holding the patterning means, the substrate 
table, the radiation system, the patterning means and the projection system, 
25 to optimize a merit function which weighs and sums lithographic errors causing anomaly in a projected image 

in the radiation sensitive layer; and 

a processor for calculating at least one compensation to be applied by said compensation means on the basis 
of a plurality of coefficients which quantify the relationship between at least one lithographic error and properties 
of at least one of the patterning means, the projection system, the radiation sensitive layer on the substrate 
30 and the projection beam. 

17. Apparatus according to claim 16, wherein said properties are selected from the group consisting of: pattern fea- 
tures: mask pattern errors; the aberration of the projection system; data concerning the radiation sensitive layer 
on the substrate; radiation beam energy and wavelength. 

35 

18. Apparatus according to claim 16 or 17, wherein said at least one lithographic error is selected from the group 
consisting of: feature displacement in the X,Y-plane; feature rotation; asymmetric magnification of feature; distor- 
tion: focal plane deformation; asymmetric distribution of CD over a target portion; asymmetry of CD with respect 
to defocus; Left-Right asymmetries of CD; L1 -L2; L1 -L5; differences of CD between patterns that are substantially 

40 directed along two mutually orthogonal directions; variation of CD within a feature, along a bar, commonly known 

as C-D; and combinations thereof. 

19. Apparatus according to claim 16, 17 or 18, wherein said coefficients comprise linear constants of proportionality 
between said properties and said lithographic error. 

45 

20. Apparatus according to any one of claims 16 to 19, wherein said coefficients comprise a plurality of families of 
coefficients for at least one of: different illumination settings of the radiation system and/or projection system; and 
different features of the pattern to be projected 

so 21 . Apparatus according to any one of claims 1 6 to 20, further comprising a database of said coefficients and/or said 
properties. 

22. Apparatus according to any one of claims 1 6 to 21 , wherein said properties comprise a field map of aberration as 
a function of position in the image to be formed at said target portion. 

55 

23. Apparatus according to any one of claims 1 6 to 22, wherein said processor is arranged to perform said calculation 
to optimize the merit function for the most critical feature of the pattern to be projected. 
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24. Apparatus according to any one of claims 1 6 to 23, wherein the compensation comprises at least one of: translation 
of the patterning means; rotation of the patterning means; translation of the substrate; rotation of the substrate; 
movement or deformation of optical elements : preferably optical elements of the projection system; change of 
energy of the radiation beam; change of wavelength of the radiation beam; change of index of refraction of gas- 

5 filled spaces traversed by the projection beam; change of the spatial distribution of the intensity of the radiation; 

change of used portion of field. 

25. Apparatus according to any one of claims 16 to 24, wherein: 

10 said properties comprise a field map of aberration of the projection system expressed in terms of the Zernike 

expansion; 

said lithographic error comprises at least one of: the X displacement of the pattern; the Y displacement of the 
pattern; the Z displacement of features of the pattern oriented in the X direction; and the Z displacement of 
features of the pattern oriented in the Y direction, each as a function of position in the X, Y plane of the target 
is portion of the substrate; and 

said coefficients are the constants of proportionality giving the sensitivity of the lithographic error to each 
Zernike aberration component in said expansion. 

26. Apparatus according to any one of claims 16 to 25, wherein: 

20 

said merit function describes the left-right asymmetry of the pattern; 

said coefficients relate the left-right asymmetry to the position of the pattern with respect to the full field of the 
projection system; and 

said compensation is adjustment of the position of the pattern with respect to the full field of the projection 
25 system. 

27. Apparatus according claim 26, wherein said compensation comprises a shift of the pattern in the X direction, such 
that the patterned beam traverses the projection system in a decente red manner. 

30 28. Apparatus according to any one of claims 1 6 to 27, wherein: 

said merit function describes the left-right asymmetry of the pattern; 

said coefficients are the constants of proportionality giving the sensitivity of the left-right asymmetry to aber- 
ration of the projection system; and 
35 said compensation is adjustment of the linear comatic aberration of the projection system. 

29. Apparatus according to any one of claims 1 6 to 28, wherein: 

said patterning means comprises a phase shift mask; 
40 said properties comprises the phase shift error of the phase shift mask; 

said merit function describes the Bossung tilt; 

said coefficients are the constants of proportionality giving the sensitivity of the Bossung tilt to the phase shift 

error and to aberration of the projection system; and 

said compensation is adjustment of the aberration of the projection system. 

45 

30. Apparatus according to claim 29, wherein said aberration of the projection system is particularly the spherical 
aberration, and more preferably the Zernike aberration component known as Z9. 

31. An apparatus according to any one of claims 16 to 30, wherein the support structure comprises a mask table for 
so holding a mask. 

32. An apparatus according to any one of claims 1 6 to 31 , wherein the radiation system comprises a radiation source. 
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Fig.2. 
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